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Measurements of angular correlations of hadrons with a (semi-)hard trigger hadron in heavy- 
ion collisions at RHIC show large angular structures opposite to the trigger which were a priori 
unexpected. These away side large angle correlations were first observed in two-particle correlations 
P, Ql and have recently also been investigated in three-particle correlation measurements 0,131 ■ We 
show that the correlation signal can be understood in terms of sonic Shockwaves ('Mach cones') 
excited by hard partons supersonically traversing the medium. The propagation of such shocks 
through the medium evolution is treated in a Monte Carlo (MC) framework :6, . We demonstrate 
that two- and especially three-particle correlations offer non-trivial insight into the medium-averaged 
speed of sound and the evolution of flow. Our findings imply that the assumption of "defiected jets" 
is not necessary to account for the observed correlations. 



PACS numbers: 25.75.-q,25.75.Gz 



The four detector collaborations at RHIC have an- 
nounced that a new state of matter distinct from ordinary 
nuclear matter has been observed in ultra-relativistic 
heavy- ion collisions The interaction of hard par- 

tons created at the first moments of the collision with 
the soft bulk medium has been suggested as a promis- 
ingprobe to study the properties of this state of matter 
[irSi E, [ill ll^, ■ Recent measurements of two- and 
three particle correlations involving a semi-hard trigger 
parton with 2.5 GeV < pr < 4.0 GeV and an associate 
threshold of 1.0 GeV < px < 2.5 GeV have shown an 
a priori surprising splitting of the away-side peak for all 
centralities but peripheral collisions which is very differ- 
ent from the broadened away-side peaks observed in p-p 
and d-Au collisions [1,]. The structure of the away-side 
peak is also persistent for different associate hadron mo- 
mentum cuts (cf. [1]). 

The observed structures have been interpreted in terms 
of interactions of the away-side parton with the soft 
medium. In particular the excitation of colorless and 
colored collective modes has been suggested, see 
[TgI . [17| . The predictions inferred from the excitation of 
sound modes [3l| turn out to be in good agreement with 
two-particle correlation data [1, Model calculations 
for Mach shocks in a true hydrodynamical framework 
have so far only been performed using a boost-invariant 
'Mach-wedge' as an approximation of the situation close 
to midrapidity ITsj. However, this is not the experi- 
mental situation in which the trigger condition only con- 
straints the near-side hadron to be at rapidity y — 0. 
The position of the away-side parton can be inferred as 
a conditional probability distribution Pjv) calculable in 
perturbative QCD. We have shown in [3] that the mea- 
sured correlation structure arises from folding of P{y) 
with the rapidity structure of the correlation for single 
away-side partons, taking detector acceptance into ac- 
count. However, as the correlation signal from single 



partons is elongated in ^, y-space by longitudinal flow 
this does not lead to a strong reduction of the predicted 
maximum correlation angle which would be in contradic- 
tion to experiment, see [7] and especially Fig. 3 therein. 

Here we focus on central collisions. In order to calculate 
two- and three-particle correlations we proceed as follows 
(for details see l6(, |7|]): We simulate the trigger condition 
as closely as possible by using a MC approach. The hard 
vertices are generated with a distribution weighted by the 
nuclear overlap and parton type and momenta are de- 
termined by randomly sampling partonic transverse mo- 
mentum spectra as generated by the VNI/BMS parton 
cascade as described in 19,]. For the energy loss cal- 
culation we use the formalism outlined in Ref. [l^] to 
obtain the probability P{AE) of the parton to lose the 
amount of energy AE from two characteristic quantities, 
namely the plasma frequency ujc and the averaged mo- 
mentum transfer {qL){YQ,(f). For the description of the 
dynamical evolution the fireball model outlined in Refs. 
pH . [23 , [23l . [2^ is used. It is constrained by transverse 
mass spectra of hadrons, Hanbury-Brown-Twiss (HBT) 
correlation radii and hadronic dN/dy distributions and 
is known to lead to a good description of direct photon 
yields as well as Raa in the high pr-region, [25, 26). It 
is also successful in predicting the back-to-back correla- 
tions for a hard trigger > 8 GeV [13] ■ Here we study 
three different cases for the dynamics of the fireball. The 
first one is referred to as 'box profile' (the model de- 
scribed in Ref. [2l|) and is characterized by a transverse 
Woods-Saxon entropy density profile with a small surface 
thickness d^s ~ 0.5 fm. This distribution is favoured by a 
simultaneous fit to transverse mass spectra and HBT cor- 
relation radii. The second utilizes the nuclear profile Ta 
as transverse density distribution as expected for a soft 
bulk matter production mechanism. In the third case 
we retain the Ta density profile but change the acceler- 
ated longitudinal expansion of [2l| into a boost-invariant 
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Bjorken evolution. The other geometrical scales of the 
model are as in the first case. Equal temperature con- 
tours for the three cases are discussed in [13] (Fig. 3 
therein). The core temperature evolution as a function 
of proper time r for central 200AGeV Au-Au collisions 
for the three different cases is shown in Fig. [T] 
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FIG. 1: Core temperature evolution in 200 AGeV Au-Au col- 
lisions from fireball models with three different transverse en- 
tropy density profiles and accelerated longitudinal or boost- 
invariant longitudinal expansion respectively (see text). 



where Cs{t) is calculated along the path of propagation. 
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FIG. 2: Calculated three-particle correlation signal as inferred 
from Mach shocks created by supersonically traveling partons 
in the medium on the away-side. A Monte Carlo simulation 
accounts for the near-side trigger condition. The transverse 
entropy density is distributed according to a Wood-Saxon 
form with 0.5 fm surface thickness ('box density'). 



We calculate the energy lost from the near-side parton 
and decide if the trigger condition is fulfilled. If the ex- 
periment triggers on a semi-hard hadron in the transi- 
tion regime where both recombination and fragmentation 
might be important, only the parton is required to ful- 
fill the trigger condition for simplicity. We have checked 
that the shape of the correlation signal as predicted by 
the model is not sensitive to the exact value of the trigger 
threshold (we see however an uncertainty in the absolute 
normalization) . If the near-side parton passes the trigger 
condition, the direction of the away-side parton is deter- 
mined taking into account intrinsic kx- In this way one 
finds for each scenario a specific distribution of initial 
vertices, directions and momenta of away-side partons, 
i.e. all trigger bias is taken into account consistently. In 
the MC, the away side correlation is inferred from a suffi- 
ciently large sample of generated events. For every single 
event the averaged away side energy loss per proper time 
dE/dr is calculated. Our main assumption is that a frac- 
tion / of the lost energy excites a collective mode of the 
medium which is characterized by a dispersion relation 
E = CsP, where the speed of sound Cs is determined lo- 
cally by the equation of state as fitted to lattice results 
[2^. This determines the initial angle of propagation of 
the shock front with the jet axis, (/> = arccoscs. Elements 
of the front are subsequently propagated with the local 
Cs and (f) is continoulsy adjusted as 
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Since the mode is propagated with Cg relative to the sur- 
rounding medium, the shock front can be significantly 
distorted in position space by local transverse and lon- 
gitudinal flow 0, d^. Note that the final longitudinal 
position Zfinai at freeze-out is relevant only insofar as 
the longitudinal flow at Zfinai determines a longitudinal 
boost in momentum space [7|. Once an element of the 
shock front reaches the freeze-out hypersurface, the sonic 
mode cannot propagate any further. We assume that 
the energy and momentum contained in the shock front 
are converted into additional kinetic energy and employ 
the matching condition to calculate the additional boost 
Mshock of a volume element received from the shock front. 
The Cooper-Frye formula is used to convert the fluid el- 
ement into a hadronic distribution: 
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As discussed in [30| section 4 the observed yield is very 
sensitive to flow if the associate hadron cut is set well 
above the bulk matter momentum scales at freeze-out 
(^ 3 T/) and comes from those parts of the shock front 
where the transverse flow is (approximately) aligned with 
the emitted particle momentum. In the presence of a 
Shockwave the signal at 1 GeV is in fact about 9 times 
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FIG. 3: Calculated three-particle correlation signal as inferred 
from Mach shocks created by supersonically traveling partons 
in the medium on the away-side. A Monte Carlo simulation 
accounts for the near-side trigger condition. Upper panel: 
transverse entropy density distributed according to the nu- 
clear profile Ta; accelatered longitudinal expansion. Lower 
panel: transverse density as above, but boost-invariant longi- 
tudinal expansion. 



distortions in position space by flow perpendicular to the 
wave propagation do not map into a sizeable angular dis- 
tortion of in momentum space since flow and shock have 
to be (approximately) aligned in every case were sub- 
stantial contributions to the yield are expected. On the 
other hand momentum conservation dictates that even if 
shock and flow are not aligned the correlation signal can- 
not vanish. In this case a broader correlation structure 
is observed at a lower pT ■ 

This is important for the understanding of 2- and 3- 
particle correlation signals: If the direction and initial 



position of the away-side parton in a given MC event is 
just such that transverse flow and cone are aligned, the 
expected correlation signal is very different from the case 
where transverse flow and the cone are orthogonal to each 
other (cf. left panel of Fig. 2 in [13] )■ In the aligned case 
both wings of the Mach cone leave their trace in the cor- 
relation signal while in the case were transverse flow is 
orthogonal to one wing of the shock cone the correlation 
signal can be suppressed and particles are redistributed 
to lower associate pt- Note that the case were only one 
wing appears in the correlation signal of a specific away- 
side jet could be misinterpreted as a "deflected jet" . The 
calculated two-particle correlation signal which is com- 
pared with the experimental result emerges in the Monte 
Carlo approach as an average of the signal expected from 
individual events. Large event by event fluctuations of 
signal shape and strength are observed. The averaging 
reduces information: if one would study a subsample of 
MC events in which always only one wing of the cone 
is visible, after averaging the signal would be difficult to 
distinguish from a situation in which always both sides of 
the cone are seen with half strength. It is therefore clear 
that an averaged two-particle correlation measurement 
can in principle not distinguish the two situations. 

However, more information is carried by 3-particle cor- 
relations. Both the PHENIX and STAR coUaboration 
have presented (preliminary) 3-particle correlation mea- 
surements d, We chose the way of presenting the 
correlation signal in the 0i2 and (pi^ plane as utilized 
by the STAR collaboration. Here 0i2 indicates the an- 
gle between the trigger hadron and one of the associated 
hadrons and (pi^ the angle between trigger and other as- 
sociated hadron. 

Our Monte Carlo approach to calculate the two-particle 
correlation signal can be straightforwardly extended to 
calculations of the three-particle correlation signal in the 
away-side hemisphere assuming that 3-particle correla- 
tions can be calculated from factorized 2-particle corre- 
lations (as we do not simulate the complete near side jet 
but only the leading hadron, we do not show correlation 
signals inside the trigger hemisphere). Within this as- 
sumption, we neglect true 3-particle correlations. Since 
the Mach cone is a collective phenomenon it can involve 
typically of the order of 20 — 50 particles, depending on 
the specific event under consideration Q. While all of 
them are strongly correlated with the original away-side 
parton (and hence with the trigger) and thus correlated 
among each other because of the conic geometry, these ef- 
fects are part of the factorized two-particle correlations. 
Additional correlations of particles inside the cone are 
weak and are only induced by momentum conservation 
T.Pl= Paway and ^ px = (see also Q). 
We show the results for the three-particle correlation pat- 
tern on the away side in Fig. 2 and Fig. 3 for the three 
different scenarios discussed above. Note that we assume 
for this calculation that the correlation measurement is 
performed exactly at midrapidity. All cases clearly indi- 
cate a correlation signal along the diagonal in the <f>i2 and 
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(/>i3 plane and variable strength along the off-diagonal. 
The maxima along the diagonal mirror the maxima in 
the two-particle correlation measurements. The opening 
angle of the Mach cone (i.e. the distance between the 
maxima) is larger whenever the averaged Cg is reduced. 
Since Cg is lowest in the vicinity of the phase transition 
[2^ . the opening angle is the largest for the Ta entropy 
density profile in the Bjorken evolution case, a bit less in 
the box entropy density case and least in the Ta entropy 
density case with accelerated longitudinal expansion, cf. 
the T evolution in Fig. 1. 

Correlation strength along the diagonal is expected if two 
correlated particles on the away-side are typically picked 
up from the same wing of the cone, while the off-diagonal 
structures appear if they originate from opposite wings. 
The three particle correlation signal for the box density 
profile has most strength on the diagonal, indicating that 
events in which only one side of the cone is aligned with 
the transverse flow at the freeze-out hypersurface occur 
more frequently. The off-diagonal correlation signal is 
more pronounced in the Ta entropy density profiles. This 
demonstrates that while off-diagonal correlation signals 
are a strong indicator for the existence of the formation of 



Mach cones on the away-side, their strength is sensitive 
to the details of the underlying evolution. We emphasize 
that it is not necessary to invoke the notion of a "de- 
flected jet" as an explanation for the diagonal correlation 
structures. 

This analysis demonstrates that the tomographic capa- 
bilities of correlation measurements are significantly in- 
creased by three particle correlations [1, 0|. They open 
the possibility to test properties of the excitation of col- 
lective modes [13, [H, [ij- Furthermore they may also 
provide a deeper insight into the speed of sound in the 
nuclear medium as well as into our theoretical under- 
standing of the evolution of the flow field and the bulk 
matter density in ultrarelativistic heavy-ion collisions. 
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